Peak shifts due to B^ — B^> rescattering in T(5S f ) dipion transitions 
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We study the energy distributions of dipion transitions T(55) to T(lS,2S,3S)n + n~ in the fi- 
nal state rescattering model. Since the T(55) is well above the open bottom thresholds, the di- 
pion transitions are expected to mainly proceed through the real processes Y(55) — > B"B'*' 
and B M B M -» T(15, 25, 3S)7T + 7r _ . We find that the energy distributions of T(l S, 25, 3S)tv+tv' 
markedly differ from that of T(55) — > B^B^ . In particular, the resonance peak will be pushed up 
by about 7-20 MeV for these dipion transitions relative to the main hadronic decay modes. These 
predictions can be used to test the final state rescattering mechanism in hadronic transitions for 
heavy quarkonia above the open flavor thresholds. 
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I. INTRODUCTION 



Hadronic transitions between heavy quarkonia are im- 
portant for understanding both the heavy quarkonium 
dynamics and the formation of light hadrons (for recent 
reviews see, e.g. [l|). Particularly, in recent years, the di- 
pion and single pion transitions have proved to be a very 
efficient way to discover new or missing charmonium or 
charmonium-like states, the so-called "X, Y, Z" mesons, 
by Belle, Babar, and CLEO collaborations (see [2| for a 
recent review and related references). To search for the 
partners of X and Y in the 66 sector ,say Xb and Y,, a 
similar approach was also suggested More recently, 
Belle has found a very striking result [J] that signals of 
T(m5)7r + 7r~ with m = 1,2,3 collected at the energy 
of 10870 MeV, i.e. at the peak of T(55), indicate that 
the partial widths of Y(55) — * T(m5)7r + 7r~ are larger 
than that of T(45) -> T(15, 2S)tt+tt- by two orders of 
magnitude or more if the T(55) is the sole source of the 
observed dipion transition events. 

If the peak at T(10870) is identified with the T(55) 
resonance, we will have to answer the question: what is 
the essential difference between the T(55) and T(45) 
in their dipion transitions. To answer this question, 
in an earlier paper [5| we use the final state rescatter- 
ing model [|| to study the dipion transitions of T(55) 
and T(45). In this model, the Y(55/45) first decays to 
B(*'B(*', and then the B meson pair turns into a lower 
T state and two pions through exchange of another B^ 
meson. Since the real rescattering contributions are ex- 
pected to be dominant in these processes, the difference 
between T(55) and T(45) dipion transitions can be ex- 
plained mainly by the difference in available phase space 
for the T(55/45) -> flWJW decays. Although sharing 
the same quantum numbers and having similar lcptonic 
widths, the Y(55) and T(45) are dramatically different 
in their hadronic transitions. In the real rescattering 
process where the 66 resonance is above the open bot- 
tom threshold, the amplitude is basically proportional to 
the probability of the corresponding open bottom decay, 



which is proportional to the P-wave phase space factor 
|Pi| 3 [H 01 • Here, p\ denotes the 3- momentum of B^ 
or i?M i n the rest frame of T (55/45) ._ Note that the 
p- value \pi\ of the decay T(55) — > BB is about 3.84 
times larger than that of T(45) — > BB, and this fact will 
mainly result in a huge difference, which is about a factor 
of 200-600 in magnitude pS}, between the partial widths 
of dipion transitions of T(55) and T(45). (Note that the 
|pi| 6 factor would give an enhancement of about 3200 for 
the T(55) relative to T(45) but an effective form factor 
associated with the coupling constants will lower the en- 
hancement factor by about an order of magnitude, see 
discussions in next sections.) As a result of the rescat- 
tering mechanism, it might be unnecessary to introduce 
an exotic interpretation of T(55) resonance or a Y, state 
to account for the experimental data 0]. 



To further clarify the issue mentioned above, it is use- 
ful to study other features of the final state rescattering 
mechanism. In this paper we will show that a distinct 
and important consequence of the final state rescatter- 
ing mechanism is the peak shift effect, which is related 
to the unique energy dependence of the cross section of 
e + e~ — ► T(55) — * Y(m5)7r + 7r~. As a result of the final 
state rescattering process B^B^ — ► T(m5)7r + 7r~, the 
energy distributions of the T(55) — > T(m5)7r + 7r~ events 
will differ from that of T(55) -► BWfiW. In particular, 
the observed resonance peak in T(55) — > T(rn5)7r + 7r _ 
is expected to be pushed up to higher energies markedly 
compared with that of T(55) -> SWfiW. The physical 
reason for this peak shift is quite evident that at higher 
energies around the T(55) resonance peak the rescatter- 
ing process will acquire more phase space and then get a 
larger rate. This peak shift effect is an inevitable result 
of the final state rescattering, and therefore it can serve 
as a crucial test for this mechanism. In the following, we 
will study this effect in a more quantitative way. 
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FIG. 1: Typical diagrams for T(nS) -> B W B (,) -> T(raS)S. 
Other diagrams can be obtained by charge conjugation trans- 
formation B <-> B and isospin transformation B° <-> B + and 
S°<-> fi- 



ll. THE RESCATTERING MODEL 



As in Ref. [H, we assume that in the T(45, 55) dipion 
transitions the two pions are produced mainly via scalar 
resonances coupled to intermediate B^ mesons due to 
the long-distance final state interactions. The typical 
rescattering diagrams for T(45, 55) — > T(15, 25)5 are 
shown in Fig. [TJ and the others can be related to those 
in Fig. Q] by charge conjugation transformation B <-> B 
and isospin transformation B° «-> _B + and S <-> £? _ . 
Here, 5 denotes scalar resonance a or /q(980) (perhaps 
also /o(1370)), which will decay to irir(KK) eventually. 

In Fig. [TJ the intermediate states B^B^ can be real 
or virtual, which corresponds respectively to the imagi- 
nary part or the real part of the amplitude. In Ref. [5j, we 
argued that in general it is the real rescattering process 
that is dominant unless the resonance is very close to the 
open flavor threshold. More quantitatively, we used the 
dispersion relation to estimate the virtual rescattering 
contributions and found that they are extremely small 
for T(55) because T(55) is far above the BB thresh- 
old, whereas for T(45), they can be comparable to the 
real ones but with large uncertainties ||. Nevertheless, 
this does not affect the calculated large difference be- 
tween the transition widths of T(55) and T(45). (In 
contrast, in the case of A(3872) jg| or Z(4430) [g, the 
virtual rescattering effects could be dominant due to the 
extreme closeness between the resonance mass and the 
open flavor threshold.) 

In the real rescattering processes, the amplitudes 
of Fig. QJa,b,c,d) are dominated by their absorp- 
tive(imaginary) parts, which can be derived by the 



Cutkosky rule as [|[ 
\pi\ 



Abs, ; = 



dtoAi{T(nS) bMb'*') 
xCi(flW5W _» T(m5)5), (1) 



where i = (a, b, c, d), and d£l and p\ denote the solid an- 
gle of the on-shell system and the 3-momentum 
of the on-shell B^*' meson in the rest frame of T(n5), 
respectively. 

The amplitudes At and Ci are determined by the effec- 
tive Lagrangians @ 
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where d = d — d . Here, the heavy quark symmetry 
and chiral symmetry in ^ are basically ensured by the 
kinematic conditions. 

Following Ref. [5j , we choose the on-shell coupling con- 
stants 

9r( m s)BMBM = 24, m < 4 (3) 
5t(5S)bs = 2.5, (4) 
9r(5S)B'B = 1-4 ±0.3, (5) 

<?T(5S)B*B* - 2.5 ±0.4, (6) 

where the values in ([3]) and ([4J6]) are determined by 
the measured widths T(T(45) -> BB) and T(T(55) -> 
£j(*)_B(*)) 5 respectively. On the other hand, the coupling 
constants gsB(*)B(*i can be related to the well-known one 
9d*Dtt [H by heavy quark flavor symmetry and chiral 
symmetry, and we choose [|[ 

g<jBB = 9aB'B' = 10 GeV, 
gjoBB = g/oB'B' = 10v2 GeV. 

To account for the off-shell effect of the exchanged i?W 
meson in Fig. [TJ one need introduce form factors, such 
as 



Fi(mi,q 



(A + rrii) 2 — m\ 

(A±TOi) 5 



r 



(7) 



to the vertexes SBB and T(mS)BB. We will fix the 
cutoff A = 660 MeV @ in our numerical analysis in the 
next section. 

We treat the scalar resonance 5 as a narrow one and 
use the Breit-Wigner distribution 
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VtT s (t) 
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m|r s (i)2 



(8) 
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TABLE I: Resonance parameters of a and /o(980) [Tol. UA 



ms 


QStvtt 




gSKK FsKK 


(MeV) 


(GeV) 


(MeV) 


(GeV) (MeV) 


g 526 ± 30 


3.06 


302 ± 10 




/o(980) 980 ±10 


1.77 


61 ± 1 


2.70 12 ± 1 



to describe the resonance in the calculation of cross sec- 
tions, as the treatment of p resonance in Ref. p|. In ©, 
the variable t denotes the momentum squared of S, and 
the function Ts{t) is given by 



r 5 (t) 

Pit = 



Ptt9Stttt . PKQSKK 



8nt 



8lTt 



(9) 



Pk 



K 



The resonance parameters in (|8|) and the coupling con- 
stants in Q, which are listed in Tab.JH are chosen mostly 
from Ref. [n|, except for to /o ( 9 80) from PDG2006 pjj . 

Although the Breit-Wigner description © for the 
scalar resonance is somehow rough, especially for the 
(7, it is efficient enough to get the order of magnitude 
of the widths r(T(45/55) -> T(mS>7r) correctly @. 
Moreover, since the final state phase space (FSPS) of 
T(55) — > T(mS)S is smeared by the large width of S, 
one would not expect that this FSPS could bring on 
large energy dependence of the cross section, unless it 
happens to be near the threshold, such as the case of 
T(55) — ► T(35)(T. We will discuss this further in the 
next section. 



III. PEAK SHIFTS IN T(55) DIPION 
TRANSITIONS 

In the absorptive part Abs^ in (JTJ) , the amplitude 
Ai is proportional to |pi| since it involves an on-shell 
P-wave vertex T(nS)B^ B^*\ Furthermore, a hidden 
factor \pi\ will emerge after performing the integral in 
(P) explicitly. As a result, the amplitude Abs^ is pro- 
portional to |p*i| 3 . As we have mentioned above, this 
brings on the strong s-dependence of the cross section of 
e+e- -> T(55) -> -> T {mS)TT+ n~ , and may 

change the distribution of the T(m>S)7r + 7r~ signals sig- 
nificantly around the peak of T(5S) resonance. 

In general, neglecting radiative corrections and the 
beam-energy spread, the cross section for the process 
e + e~ — > resonance r — > hadronic final state / at 
the center-of-mass energy y/s can be approximately ex- 
pressed by the Breit-Wigner form 



a(s) = 12tt 



r--r/( s ) 



(s — 77T, 2 ,) 2 + m 2 r r (s) 2 



(10) 
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FIG. 2: Resonance line-shapes of the naive Breit-Wigner dis- 
tribution (dashed line) and the one with Y s r (s) 
BB as the intermediate state (solid line) 
are normalized to one, respectively. 



n s 



and 

The maximums 



where the resonance r is parameterized by its mass m r , 
the total width r r , the electronic width T^, e , and the par- 
tial width r£ for the decay channel /. For the T(5S), in 
(TIT)]) we have neglected the weak s-dependence of m r , T^ e . 
However, for the final states / = Y(lS/2S/3S)ir + Tr~ , the 
s-dependence of T(. is as strong as |pi(s)| 6 , and therefore 
can not be neglected at all. 

As a first step, we will neglect the energy dependence of 
the total width of T(55) as well. Thus, for the final state, 
such as / = light hadrons, the energy distribution of the 
cross section of T(55) will peak at m^s exactly as that 
shown in Fig.[5]by the dashed line. We call this the naive 
Breit-Wigner distribution. In Fig. [2J the resonance pa- 
rameters of T(55) have been chosen from PDG2006 [l(| 
with the following central values of m 5 s and r 5 s: 



m 5S = 10865 MeV, 



T 5S = 110 MeV. 



(11) 



Then, we can compare the naive Breit-Wigner distribu- 
tion with the one where T f 5S (s) ~ |pi(s)| 6 . For the BB 
intermediate state, the result is shown in Fig. [5] with the 
solid line. One can see that the strong s-dependence of 
r^g pushes the resonance peak up with a energy shift of 
about 15 MeV. 

The two-body decay modes B^B^ and B^ B S W , 
which are dominant ones of T(5S) [10|, will bring on 
large energy dependence of T§s and change the distri- 
butions in Fig. [2] consequently. For simplicity, we will 
choose the main decay mode B*B* to estimate the en- 
ergy dependence of T$s- That is, 



r 5S (s) 



5S 



3m|» 



\/s — 4m 2 B , \ 3 



(12) 



which can be derived from the effective Lagrangian (|2c|) 
with .9t(5S)s*b* being treated as a constant independent 
of s. Here, T 5 s(mss) is normalized to be 110 MeV, like 
that in (flTj). 
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FIG. 3: The energy dependence and a dependence of T5s(s). 
The dashed line, solid line and dotted line are evaluated with 
a = 0.0, 0.6 and 1.0 GeV -2 , respectively. 



FIG. 4: Resonance line-shapes for / = B* B* (dashed line) 
and the one for / = T(lS)7r + 7r _ (solid line) and for / = 
T(3S)tt + tt~ (dotted line) with a = 0.6 GeV" 2 . The maxi- 
mums are normalized to one, respectively. 



Furthermore, one need introduce form factors for every 
T(55)BW J BW vertexesinFig.Hl Following Ref. [3, we 
choose the form factor 



Exp(-a|pi(s)| 2 ) 
Exp(-a|pi(™ 2 s )| 2 ) 



(13) 



which is normalized to 1 at s 



'5Si 



since we have ex- 



tracted the on-shell coupling constants 3t(5S)b(*)b(») m 
(jUini) at this energy point. This form factor, which is as- 
sociated with the effective T(5S)B^ coupling, can 
be understood as the consequence of the overlap integral 
of the wave functions of T(55) and the two B^ mesons, 
and the scale factor a can be related to the effective ra- 
dius of the interaction, R, by a = R 2 /6 [12] • It is this 
form factor (or others of a similar form) that plays the 
role to naturally balance the otherwise over-increased de- 
cay rates with increased phase space. Namely, this form 
factor is a sort of " cutoff' for the high momentum of 
the final state mesons. The momentum dependence of 
the effective couplings can partially explain why the cou- 
pling constants 3t(55)b(*)s<*) are smaller than the one 
9r(4S)BBt though for a = 0-1 GeV~ 2 the suppression 
from the form factor is not strong enough to account for 
the large difference between the values in © and in (0J 
[6]). This suppression may be additionally due to the node 
structure of the wave function of highly excited T(55) [5J. 
Needless to say, the form factor in l|13p should also 

change the energy dependence of the width Tfg ' (s) 
and the total width rss(s) in (fT2)) subsequently. Such a 
dependence of Fss(s) is shown in Fig. [31 One can see that 
the energy dependence is significantly weakened when a 
is about 0.6 GeV -2 (solid line in Fig. [3J. As a result, 
the distribution of T(5S) with / = B*B* and a = 0.6 
GeV -2 , which is shown in Fig. H] with the dashed line, is 
very close to the naive Breit-Wigner one. 

With a = 0.6 GeV" 2 QJ], the distributions for / = 
T(lS)n + Tr~ can be evaluated numerically. Since here we 



only focus on the line-shape and the peak shift, which 
are almost independent of the parameters except a and 
r5s('7i5s), it is quite safe to choose the central values 
for other parameters introduced in the last section. The 
distribution for / = T(1S)tt + tt~ is illustrated in Fig. [4] 
with the solid line. Compared with the distribution for 
/ = B*B* (the dashed line), the peak shift is about 7 
MeV. If we choose the channel BB as the main hadronic 
channel, the peak shift can be as large as 11 MeV for the 
T(1S)tt + it~ channel. 

The energy distribution for / = Y(2S)7r + 7r~ has sim- 
ilar line-shape to that for / = T(1S)tt + it~ within the 
energy region of 10780-10950 MeV. Furthermore, since 
the transition T(55) — ► T(3S)tt + tt~ is also observed [3| 
with a large rate and quite high statistical significance 
(3.2a), we also evaluate its energy distribution in Fig. [4] 
with the dotted line. We find there is an additional 
peak shift of about 12 MeV relative to the distribution 
for / = Y(15)7r + 7r~. This is just because the role of 
the scalar resonance in Fig. [T] In our model ||, the 
final-state 7r + 7r~ are assumed to be dominated by the 
scalar resonance S(cr, /o(980)...) [Tlj | . which can be de- 
scribed by the Breit-Wigner distribution {8} at the cross- 
section level. If the mass of the a resonance is chosen 
as m a — 526 MeV given in Tab. HI then one can eas- 
ily find that the mass difference m^s — m 3S ~ 510 MeV 
happens to be in the center of the distribution in {8]). 
As we have mentioned in the end of the last section, 
this brings on another strong s-dependence to the cross 
section for / = T(3S')cr(7r + 7r _ ). Thus, the distribution 
for / = T(3S)tt + it~ is pushed up farther than those for 
/ = T(15, 2S)ir + ir~. Similarly, there will be a long tail 
in the distribution for / = T(25)7r + 7r _ within the energy 
region of 10950-11100 MeV due to the emergence of the 
scalar resonance / (980). 

It is worth emphasizing that the scalar resonance dom- 
inance is just a simplification to study the effects of 
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T(55) in the T(15)7r + 7r _ channel compared with that in 
the B*B* channel. Similar shift should be obtained in 
the Y(25)7r + 7r~ channel. On the other hand, the shift 
in the T(35)7r + 7r~~ channel could be 12 MeV larger than 
those in the T(15, 2S)ir + ir~ channels, but this depends 
on the assumption that the scalar resonance a dominates 
the Y(35)7r + 7r~ transition. Furthermore, if the BB is 
chosen as the main hadronic channel, the peak shifts will 
be enlarged by 3-4 MeV. 



FIG. 5: a-dependence of the peak shifts compared with 
m 5 s = 10865 MeV for / = B* B* (dashed line), / = 
T(lS> + 7r- (solid line) and / = T(3S)7r + 7r _ (dotted line). 

the final state rescattering mechanism in the Y(55) — > 
T(1S,2S,3S)tt + tt" transitions [5J. Thus, the differences 
between the distributions for / = Y(15, 25, 3S)7r + 7r~ 
tend to disappear if the non-resonance contributions and 
their interference with the scalar resonances are impor- 
tant. Therefore, to detect the differences between the dis- 
tributions for / = T(15, 25, 3S)ir + Tr~ can provide useful 
information on the role of the scalar resonances in these 
transitions. More valuable information can also be given 
by the measurements on the Minn) spectrum and the 
distribution of the helicity angle [4{. 

It is obvious that the position of the resonance peaks 
of Y(55) in Y(15, 2S, 3S)ir + ir~ channels depend on the 
value of the scale factor a in the form factor (|13p . How- 
ever, the energy distributions of Y(55) in other main 
hadronic decay channels (e.g., BB...) also depend on 
the same factor a. So, it is essential to compare the line- 
shape for e.g. / = Y(15)7r + 7r~ with those main hadronic 
decay channels with different choices of the value of the 
scale factor. For a = 0.0-1.0 GeV~ 2 , the results are 
shown in Fig.0 where the B*B* channel is used to serve 
as one of the main hadronic channels again. 

In Fig. [51 with a reasonable choice for the cutoff factor 
a, say, 0.3-0.8 GeV~ 2 , the relative peak shifts for / = 
Y(15, 35)7r + 7r~ are found to be almost independent of 
the cutoff. Thus, we incline to conclude that there should 
be a peak shift of about 7 MeV for the distribution of 



IV. SUMMARY 

In summary, we study the energy distributions of dip- 
ion transitions Y(55) to Y(15, 2S, 3S)ir + ir~ in the final 
state rescattering model. Since the Y(55) is well above 
the open bottom thresholds, the dipion transitions can 
proceed through the real process Y(55) — > B^B^*' and 
the subsequent process £(*)£(*) -> Y(15, 25, 35)7r + 7r-. 
This model can not only explain the observed unusually 
large rates of T(55) -> T(15, 2S, 35)7r+7r-0, but also 
predict a unique energy dependence of the cross sections 
of e+e- -> T(55) -> T(15, 25, 3S> + tt- . We find that 
the energy distributions of T(15, 25, 35)7r + 7r~ markedly 
differ from that of T(55) -> £(*)£?(*), and in particular, 
the resonance peak will be pushed up by about 7-20 MeV 
for these dipion transitions relative to the main hadronic 
decay modes. These predictions can be used to test the 
final state rescattering mechanism in hadronic transitions 
for heavy quarkonia above the open flavor thresholds. 
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